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I. INTRODUCTION 



In a recent paper, quantum mechanical systems with explicit time dependence were 
studied via separation of variables 0. In that paper, the authors identified changes of 
variables that allowed certain time-dependent quantum mechanical systems to be treated 
using separation of variables. (For convenience, we will use the term separable to describe 
the situation in which separation of variables can be applied effectively.) The applicability 
of separation of variables in these models means that for these time-dependent systems, 
there exists a mathematical structure and formalism analogous to that of the usual time- 
independent Schrodinger equation. As an application of their approach, the authors used 
separation of variables, in combination with both operator and power series methods, to 
solve exactly the harmonic oscillator with inverse frequency linear in time. 

That paper left open some important questions. First, can the technique of be used 
to identify exactly soluble time-dependent generalizations of other exactly soluble time- 
independent systems? Second, how special was the change of variables that paper used 
to find separable models? Would other changes of variables have yielded other models? 
Third, can we use this method to simplify or make possible other exact calculations in 
time-dependent quantum mechanical theories? And finally, what connections can be drawn 
between the use of separation of variables and other approaches to time-dependent quantum 
mechanical systems? This paper addresses these questions. 

Our central findings are as follows. First, we show explicitly how to construct exactly 
soluble time-dependent generalizations of any exactly soluble time-independent model. We 
do this using a time-dependent redefinition of the spatial coordinate which is linear in the 
spatial variable. Second, we show that this linear change of variables, although apparently 
quite special, is actually quite general: if we consider any time-dependent re-definition of 
the spatial variable (but leave the time coordinate as the second variable), exactly the 
same models are uncovered as by the simple linear transformation. Including arbitrary 
multiplicative transformations of the wavef unctions, we are able to obtain a larger family of 
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separable and exactly soluble models. We find that separability singles out a unique possible 
wavefunction transformation, and makes possible a larger (but still specified) set of useful 
changes of variables. Third, as an application of our method, we present exact calculations 
of propagators in exactly soluble time-dependent models, combining our techniques with 
those of shape invariance. And, finally, at various points in the paper, we comment on 
the connections between our work and the existing literature on time-dependent systems, 
with a focus in one section on establishing the quantum generalization of the work of f|| 
on classical Hamiltonians with an invariant quadratic in momentum, and showing that this 
leads to exactly the same set of models as separability does. 

One of the key results of this paper is that we reduce the determination of whether a 
time-dependent potential is exactly soluble to the determination of whether an associated 
time-independent potential is exactly soluble. We do this through the use of separation of 
variables as the central aspect of our approach. As is well-known, this method reduces a 
multi-dimensional problem to two or more lower dimensional ones. The method of separation 
of variables is very powerful, and so has been applied broadly. As an excellent recent 
example, Sklyanin has introduced the so-called functional Bethe ansatz as an alternative 
to the algebraic Bethe ansatz, in order to remove the restrictions inherent in the second 
approach. This new method is a hybrid of the quantum inverse scattering method and the 
method of separation of variables ||. 

The existing literature on time-dependent quantum mechanics already covers a variety of 
other approaches. A sampling of such articles includes , as well as the references therein. 
These articles provide a variety of approaches to time- dependent quantum mechanics, for the 
most part with a focus on exact solubility from the start. In this way our work is somewhat 
different; we take separation of variables as our central criterion, and study its implications. 
Our interest in separation of variables is twofold: separability implies the existence of a 
mathematical structure in a theory analogous to the structure associated with the ordinary 
time-independent Schrodinger equations; and, as we see from our work, particularly as 
it stands in relation to the work in f§-|7| (e.g., our methods identify all the exactly soluble 
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models found with other approaches in those references), separation of variables also appears 
to be an appropriate stepping stone in at least one path toward uncovering and classifying 
the exactly soluble time-dependent models. The agreement among these various papers 
(including ours) where they overlap suggests the existence of some fundamental principles 
that connect their various approaches. We will comment on the connections between our 
work and that of these other papers where relevant. 

For the sake of simplicity, we will restrict our attention in this paper to systems with only 
one spatial dimension, but it will be readily apparent that our techniques can be applied 
without modification to higher dimensional models as well (see also 0). 

Overall, we feel time-dependent quantum mechanics is still inadequately understood. For 
example, an explicitly time- dependent property analogous to shape invariance — which has 
greatly enriched our understanding of exact solubility in the time-independent case — is 
still elusive. We believe the restrictive nature of the results we have obtained using changes 
of variables with separation of variables is significant, because this helps us to determine 
what the actual range of exactly soluble time-dependent models is; this can serve to focus 
the search for a fundamental explanation of exact solubility in the time- dependent case. 



II. EXACTLY SOLUBLE TIME-DEPENDENT POTENTIALS 

To begin, we consider non-relativistic quantum mechanics in the Schrodinger picture. 
The Schrodinger equation in position space reads 

H*0M) = , (2-1) 

where the Hamiltonian, with time- dependent potential, is 

Our goal is to identify time-dependent potentials V(x, t) such that the Schrodinger equation 
(27T) is exactly soluble. Obviously, separation of variables in terms of x and t will not be 
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helpful in this case, so we attempt to find a change of variables for which separation of 
variables will be useful. To this end, we define a new variable y = f(t)x. In order that this 
change of variables be productive, we require that the potential factorizes in terms of these 
new coordinates, that is, we require that 

V(x(y,t),t) = g(t)V(y) + h(t) . (2.3) 

The time-dependent Schrodinger equation now reads 

~^ f2{t) w^ {y ' t] + ( 9{t) V{y) + h(jt) ) * {Vi t] = lh W) y h^ {y ' t] + ih ^^ ■ {2A) 

We use \& to denote the wavefunction as a function of y and t; thus t) = ^/(y(x,t),t). 
Note that in (|2.4j) , the partial derivatives with respect to time are now taken with y fixed 
rather than with x fixed. 

Inspecting ( [2.4|) , it is easy to see that separation of variables in y and t will work provided 
that 

f(t) = (pt + q)- l/2 , (2.5) 

and that 

g(t)=af(t) = aipt + q)- 1 , (2.6) 

where a, p, and q are constants. The constant a may be absorbed in V(y); thus, without loss 
of generality, we take a = 1. Hereon, we take the separation of variables conditions ( ^.3[ ), 
(|2.5D, and ( P7B| ) to hold. Now by considering solutions to the Schrodinger equation ( |2.4| ) of 
the form ^/(y) = T(t)Y(y), we see that the general solution to the Schrodinger equation can 
be written as 

%,t) = ]T c k T k (t)Y k (y) , (2.7) 

k 

where the c k are arbitrary constants, and where Y k and T k are the solutions of, respectively, 

h 2 d 2 Y k Hp dY k ~ 

-2^W^ y ^y~ V{y)Yk = lkYk m 
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and 

dT 

i% (pt + q ) _± _ (pt + g )/i(t)T fe = 7fc T fc , (2.9) 

with 7/c a complex constant. From ( |2.8| ), we see that Im^) = —(hp/A). The differential 
equation for T k can be solved exactly, 

T k {t) = e-^^)( q + pty {ilk/hp) , (2.10) 

where 

rj(t) = f dth(t) . (2.11) 



o 

Thus, for the above time-dependent Hamiltonians, one recovers a mathematical structure 
formally analogous to the one one finds for time-independent systems. Defining the (non- 
Hermitian) pseudo-Hamiltonian 

~ h 2 d 2 hp d ~ , , . 

we see that the problem of solving the time-dependent Schrodinger equation is equivalent 
to finding the eigenfunctions and (complex) eigenvalues of H. (Note that in the limit that 
p = 0, the pseudo-Hamiltonian becomes the ordinary Hamiltonian, as it must.) 

We are now in a position to determine the exactly soluble potentials as follows. Defining 
the wavefunction 

Y k (y) = e- im ™ 2 / Ah Y k (y) , (2.13) 

we find that Y k satisfies an ordinary (time-independent) Schrodinger equation 

h 2 d 2 Y k . 
2m dy 2 

where 



V-(y)Y k = e k Y k , (2.14) 



V_(y) = V(y) - ^V, (2.15) 



and where is a real constant given by 



efc = Ik^ 



Tip 



(2.16) 



Thus, whenever V- (y) is an exactly soluble potential of ordinary time- independent quantum 
mechanics, then the original time- dependent system is itself exactly soluble. Hence whenever 



V(x,t) = f{t) 



VM(t)x)^f 2 (t)x* 



pt + q 



x ^p 2 m x 2 



+ h{t) 
+ Kt) , 



(2.17) 



' y/pt + q 8 pt + q_ 

and V-(y) is an exactly soluble potential for the ordinary time- independent Schrodinger 
equation, the time-dependent Hamiltonian of the Schrodinger equation ( |2.2|) is exactly sol- 
uble. 

In this way, we have identified a family of exactly soluble time-dependent potentials to 
generalize any exactly soluble time-independent potential. Of course, much work has gone 
into the study of exactly soluble time-independent potentials; most recent work has focused 
on the uses of shape invariance [§]. In fact, all the well-known exactly soluble potentials of 
quantum mechanics can be understood to derive their exact solubility ultimately from the 
property of shape invariance. If a potential V-(y) is shape invariant, this means that 

n 



V-(y) = W\y,a) - 



'2m 



W'(y,a) 



(2.18) 



and W(y, a) satisfies the condition 

n 



W 2 (y,a) 



W'(y,a) = W 2 (y,a(a)) 



h 



W'(y,a(a))R(a) 



(2.19) 



/2m \/2m 

for some functions a(a) and R(a). Using the theory of supersymmetric quantum mechanics 
0, the "levels" jk are given by 



.hp 

7o = -x T 



fc-i 



Ik 



% i'=l 



a 



3) i 



(2.20a) 
(2.20b) 
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where otj = c(a 3 -_i) and «o = a. 

One of the intriguing features of the exactly soluble time-dependent potentials we have 
identified is that, in addition to containing a term of essentially the same form as the original 
potential, they all contain an additional time-dependent harmonic oscillator term. We give 
a simple way to interpret this term below. But first let us point out that because of the 
appearance of this term, for the harmonic oscillator, and only for the harmonic oscillator, 
will the time-dependent generalization simply involve introducing a time dependence into 
the original coefficient of the potential. It is for this reason that [p]] discovered no exactly 
soluble systems other than the time-dependent harmonic oscillator. 

As promised, we give here an instructive interpretation of the harmonic oscillator term 
that appears in ( |2.17|) . Using the definition of VL (y) in ( 2.15|) , one can write the Schrodinger 
equation in y-space (|2.8| ) in the following form: 

, V 

... -11 

2m \ dy 2 



1 f . lTl ±-P^y) Y k V4y)Y k = e k Y k . (2.21) 



In other words, we can view ( |2.8| ) as describing the coupling of a quantum mechanical 
particle of mass m to a background electromagnetic field described by the two-potential 
(ip(y),A(y)) = (V-(y), (pm/2) y), and this necessitates the harmonic-oscillator-type term. 
Invoking the gauge transformation A(y) — > A(y) — (3y, one can obtain the eigenfunctions 
of the gauge-transformed Hamiltonian from those of the original Hamiltonian by observing 



that under this transformation, the eigenfunctions transform by 

Y k (y) ^ e^Q^-y)Y k (y) . (2.22) 



In this way, then, the solutions to the eigenfunction equation (|2.21 ) can be obtained from 



those of the corresponding equation with an ordinary kinetic term. We point this out because 
we suspect that proper use of such a gauge picture may yield a more fundamental way to 
understand and to classify the exactly soluble time- dependent models. 
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III. CHANGES OF VARIABLES: THE POTENTIAL 



The above set of exactly soluble models with the potentials of ( j2.17| ) arises from a con- 
sideration of the particularly simple change of variables y = f(t)x. In order to begin to 
develop a classification of separable and exactly soluble time-dependent models, it is clearly 
essential to consider other transformations, to see what separable and what exactly soluble 
models these produce. We have done this, and the answer we find is surprising. Consider 
any change of variables of the form y = P(x,t), while keeping t as the second coordinate. 
Then, neglecting Galilean transformations (as these are essentially trivial), we find that such 
a transformation can produce separation of variables only when the function P is in fact a 
function not of x and t arbitrarily, but only of the ratio xj y/pt + q, that is P(x/ y/pt + q). 
(Of course, the potential must also separate appropriately.) Furthermore, the only exactly 
soluble models uncovered by such transformations are exactly the exactly soluble models 
uncovered in ( |2.17|) by the very simplest transformation we considered. 



We now offer a proof of this result. The argument is rather technical, and the reader 
interested primarily in our results might comfortably omit this section on a first reading. 

We will build towards the general proof by first giving three preliminary results. First, 
note that making a transformation y = P(x/y/pt + q) for any function P, will, given an 
appropriate potential, lead to separation of variables in terms of y and t. Moreover, when 
the time-dependent potentials are determined which are exactly soluble using separation of 
variables and this transformation, one finds that they are exactly the same as the potentials 
already uncovered in the previous section by the linear transformation. The calculation 
proceeds exactly as in the preceding section; in the interest of brevity, we will not present 
the details here, as they are straightforward to obtain. (It is helpful in this calculation to 
invert the transformation formally via x = {\/pt + q)P~ 1 {y)-) 

Second, consider a transformation slightly more general than the linear transformation 
of the preceding section, namely the re-definition of variables 

y(x,t) = f(t)(x-u(t)) . (3.1) 



Going through the same calculation as we did for the case u(t) = 0, we find that separation 
of variables is only possible for / 7^ when, as before, 

f(t) = - ? ^= , (3.2) 
Vpt + q 

and when 



u(t) = a^Jpt + q + c . (3.3) 

(When / = 0, the only solution is a Galilean transformation; since this is a trivial possibility, 
we ignore this case.) The additional parameters a and c are easy to understand if we re- write 



( p.ip using the results for f(t) and u(t), obtaining 

y = f(t)(x-c)-a. (3.4) 

We see that the parameter a corresponds to a shift in the origin of y, while c corresponds 
to a shift in the origin of x. Thus these parameters give us no new models. Any potential 
for which separation of variables obtains when a = yields separation of variables for any 
value of a, and so allowing the parameter a adds no new models to our list ( |2.17| ). The 
parameter c also gives us no new models. If separation of variables works for V(x,t), it 
will not generally work for V(x — c,t) in terms of the same separating variables. However, 
V(x,t) and V(x — c,t) correspond to the same physical model. Including the parameter 
c thus allows us to use separation of variables directly no matter what origin one uses for 
x; the earlier transformation found all the same models, but entailed choosing a particular 
origin for x. 

The next result is a slight generalization of the preceding one ( |3.1| ). Consider now a 
transformation of the form 

y = ip(t)(x-u(t)) n . (3.5) 

Such a transformation can lead to separation of variables only when ip(t) = (pt + q)~ n ^ 2 (i.e., 
f(t) n , where f(t) is as above) and u(t) = a^pt + q + c, as before. Rather than going through 
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the same type of calculation again, one can obtain this result from the preceding one. If 
separation of variables holds in terms of y and t, then it also holds in terms of z — y l ^ n and 
t. In this way, we can reduce the problem of finding useful transformations of the type in 
( p.5| ) to finding useful transformations of the type in (|3.1| ), which we already know. 

We are now in a position to obtain our general result, namely that any transformation 
of the form y = P(x,t) leads to exactly the same separable and exactly soluble models as 
we have already obtained. To begin, we will write the function y = P(x,t) as a product, 

y = G(t)l[(x-u 3 (t)) nj . (3.6) 

For the purposes of our discussion, we take the rij to be positive but not necessarily integral, 
although the restriction to positive rij is not actually necessary. The Uj must all be distinct. 
For convenience, we define iV = Ylj=i n j- We note that any P which is polynomial in the 
variable x can be written in this form. Since our proof will hold for all polynomials of any 
order, it will also hold in the limit of an infinite but converging power series, which can 
be understood as the limit of such polynomials. The set of transformations which we can 
handle based on ( |3.6| ) is thus a superset of those functions which have a power series in x. 
In fact, the form of fl3.6p shows us that our argument covers all functions except those with 
extreme singularities. 

The essential idea behind our proof is to take the requirement that separation of variables 
hold under this transformation, and look at the implications of this requirement in different 
regions of space and time. This will enable us to place restrictions on the possible forms of 
G(t) and the Uj(t). 

First, consider this transformation at large values of x, when x » Uj(t) for all j. In 
this region, the transformation is essentially 

y « G(t)x N . (3.7) 

From our consideration of ( |3.5| ), we know that separation of variables here requires that 
G(t) oc (pt + q)~ N ^ 2 - Without loss of generality, we may re-scale the transformation so that 
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Next we consider the case that x is in the neighborhood of Uj(t). Here x — Uj(t) is a 
small parameter. Expanding to lowest non-trivial order in this parameter, we see that in 
this neighborhood 

y^G,(t)(x-u J (t)) nj , (3.9) 

where Gj(t) is 

k 

G 3 (t)=G(t) J] ( Uj (t) - Ui (t)) ni ■ (3-10) 

In this neighborhood, the transformation takes the form of ( |3.5|) , and so we see immediately 
that we must be able to write each Uj as 



Uj(t) =a j ^t + q j + c j . (3.11) 

Also, we see that Gj is 

G^) = (^=T ■ (3-12) 



Thus 



n (u^-u^y* = {-^=) nj vt+- q N . (3.i3) 

Since we know the general form of the u/s, we know that the left side of the preceding 
equation ( 3.13Q has no divergences, and thus the right side must have no divergences. This 



requires that qj = q in the expression for Gj, which in turn dictates that qj = q in each 
of the Uj, since, as we saw back in (|3.5|), the square root that appears in u(t) is the same 
square root as appears in <p(i). 

Putting all these results together, the equation (|3.13|) now reads 



k 

II ((% - a k )Vt + T + (Cj - c k )) ni = (bjpVt + T N ~ n3 ■ (3.14) 
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In order that the multiplicity of zeroes at t = — q be the same on both sides of the equation 
(and that there be no zeroes at any other value on the left side, since there are none on the 
right side), we must have Cj = Ck for all values of j and k, and aj 7^ for j 7^ k. (This is 
consistent with our requirement that all the u/s be distinct.) We define c by Cj = c for all 
j- 

We now see that in order for separation of variables to work, the transformation in (|3.6|) 
cannot be the most general product possible, but must necessarily be of the form 

1 k 

y= / n Ii( x - ajVt + q- c) nj . (3.15) 
y/t + q 3=1 

We are free to redefine the origin along the x-axis; this does not change the physical model. 
We therefore choose the origin such that c = 0. Factoring out the square roots (and recalling 
that iV = Y^j=i n j)i we see that the most general possible transformation which has the 
possibility of producing separation of variables is of the form 



Such a transformation is not a function of both x and t separately, but rather is a func- 
tion only of the combination xjyjt + q. As we discussed above, the separable and exactly 
soluble models that are identified via separation of variables following a transformation 
y = P(x/y/t + q) are exactly the separable and exactly soluble models found using our 
original linear transformation. 

Thus by considering an arbitrary coordinate transformation which leaves time as the 
second coordinate, we obtain no models that we did not already identify by considering the 
simple linear transformation. This result is quite restrictive, and it suggests that there is 
a rich mathematical structure underlying both the applicability of separation of variables 
to and the occurrence of exact solubility in time-dependent quantum mechanics. We would 
like to mention that we have tested our general proof by examining a number of specific 
changes of variables. In every case we have checked (which includes a number of polynomial 
and rational functions for which y(x) can be inverted in closed form), we have verified that 




(3.16) 
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the only changes of variables that can produce separability are the ones identified by our 
general argument. 

Thus we have shown that separation of variables only arises for transformations y = 
P(x,t) if P is a function of x/y/pt + q. Let us re-emphasize here the twofold nature of 
the consequences of this result: (1) there are no models that admit separation of variables 
using an arbitrary transformation for which a linear transformation does not suffice; (2) 
there are no exactly soluble models that we find using separation of variables and an ar- 
bitrary coordinate transform that we did not find using separation of variables and the 
linear coordinate transform. Thus, unless we consider changes of variables in which the time 
variable is redefined (a possibility we do not explore in this paper), we have, by using the 
simplest time-dependent transformation, already identified the most general exactly soluble 
time-dependent generalizations of the exactly soluble time-independent models which are 
accessible by means of separation of variables while using the original wavefunction. In 
the next section we include a multiplicative transformation on the wavefunction along with 
changes of variables, which enables us to find a somewhat larger class of separable models. 
Our results here (and in the next section) demonstrate that the space of time-dependent 
exactly soluble models is extremely constrained, and thus suggests that separability may 
be a useful tool in classifying such models and in identifying the structure underlying their 
solubility. 

IV. CHANGES OF VARIABLES: THE WAVEFUNCTION 

The fact that the series of potentials found by the method described above closely resem- 
bles the series of potentials discussed in || is very intriguing. It suggests that there should 
be a modification of the above prescription such that the potentials of @ - - and perhaps 
more — are exactly reproduced. Here we present such a modification, which also serves to 
help us ascertain the scope of models one can address via separation of variables. To this 
end, we define a new variable 
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y = f(t) (x - a(t)) 



(4.1) 



but also include now, in addition, a transformation of the wavefunction 



*(y,t) = E(y,t)e 



(4.2) 



where ^ denotes the wavefunction as a function of y and t, i.e. ^(x, t) = $f(y(x,t),t). 
(This ansatz is motivated in part by an attempt to generalize the transformation in (|2.22|) .) 
Without loss of generality we write the potential in the form 



V{x{y,t),t) = g(t)V(y) + U(y,t)g (t) . 



(4.3) 



Furthermore, since we are interested in finding new models, we only consider the case that 
V(y) is not exclusively a quadratic polynomial in y. 

Using the above definitions, the Schrodinger equation becomes 

h 2 d 2 ~dz 



2m dy 2 dy 



h 2 d§ f(t) .,a 

— — — %Ti y ~f~ %Ti — 

m dy f 3 (t) f 



ih <9$ 



h 2 d 2 $ h 2 (d$\ fa f(t) \ <9$ 



-[g(t)V(y) +U(y,t)\ 



g (t) „ ih dZ 



(4.4) 



/2(t)L *v , w n /2(t) /2(t) m 

In order that the coefficient of cfE/dy be a function only of y (which is necessary for sepa- 
rability), we must have 

,2 



$(y,t) = a(t)+b(t) V - + c{t)y + d{y) 



(4.5) 



with 



b(t) 



im f 



im a 



(4.6) 
(4.7) 



Shifts of b(t) and c(t) by a constant from these values are allowed consistent with separation 
of variables, but such constants can just be absorbed into re-definitions of a(t) and d(y). 
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Since our goal in (|4.2|) is to make separation of variables possible where it previously 
was not, we see that the choices one makes for a(t) and d(y) are arbitrary, as these will 
not change whether the model is separable. Hence, without loss of generality, we make the 
convenient choices d(y) = and 

a (t) = -ln/^ f dsa 2 (s) - % - f ds g (s) . (4.8) 



2 2k Jo w h Jo 

Thus we see that given a transformation to new variables (|4.1| ), there is an essentially 
unique $>(y,t) which has the potential to be useful in transforming the wavefunction to 
obtain separability. Using the value of Q(y,t) we have now determined, we see that the 
Schrodinger equation ( |4.4| ) simplifies to 



h 2 d 2 z. 



2m dy 2 



.2 



a m 2f 2 -ff 
1 -[g(t)V(y)+U(y,t)\ = -£rM ■ ( 4 -9) 



/ 2 (t) ww " J n P(t)dt 

Since V(y) is not simply a quadratic polynomial, in order for separation of variables in 
y and t to work, the function g(t) must satisfy 

g(t) = p(t) . (4.10) 

We now only have left to consider f(t) and a(t). Note first that if the potential factorizes 
in the sense of ( fOD , so that U(y,t) = 0, separability will only obtain provided that 

y/Xt 2 + fit + V 

and with a fixed by requiring that the coefficient of the term linear in y in ( |4.9|) be in- 
dependent of t. This case makes contact with and generalizes the results of the previous 
sections; more significantly, it reproduces exactly the set of models investigated in [Q. So 
in the special case {7 = 0, there are unique transformations (which we have found) of the 
spatial coordinate and the wavefunction which can produce separation of variables. 

In the more general case that {7^0, there is a somewhat larger family of separable (and, 
in turn, a larger family of exactly soluble) models. Separability requires that U (y, t) exactly 
cancel the y and y 2 terms, leading to 
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U{y,t) = -m—y - — — y . (4.12) 

(We can shift U(y,t) by any function of the form Uo(y)f 2 (t) or T(t) consistent with separa- 
bility, but these can just be absorbed into redefinitions of V(y) and go(t), respectively) 
Now by considering solutions to the Schrodinger equation (|4.9| ) of the form E(y,t) = 



T(t)Y(y), we see that, with U(y,t) as in ( |4.12| ), the general solution to the Schrodinger 
equation can be written as 

E(y,t) = c k T k (t)Y k (y) , (4.13) 

k 

where the c k are arbitrary constants, and where Y k and T k are the solutions of, respectively, 

-W^Y^^ (4.14, 

and 

ih^ = e k f 2 (t) . (4.15) 
The differential equation for T k can be solved exactly, 

T k (t) = e-^ ekT{t) (4.16) 

where 

r(t) = t dtf(t) . (4.17) 



o 



Thus, whenever V(y) is an exactly soluble potential of ordinary time-independent quan- 
tum mechanics, then the original time- dependent system is itself exactly soluble. Grouping 
togehter several terms by defining 

h(t) = - m \ ^ ~ oi 2 - act] + g (t) , (4.18) 



2/2 



we see that whenever 



m 2f 2 - // 
2 p 
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x 2 + h(t) , (4.19) 



the model is separable, and that when, in addition, V(y) is a standard exactly soluble 
potential for the time-independent Schrodinger equation, the time-dependent Hamiltonian 
of the Schrodinger equation ( |2.2|) is exactly soluble. In this way, we have identified a more 
general family of exactly soluble time-dependent potentials. 



V. THE LEWIS LEACH APPROACH 



The series of potentials ( [4.19|) that we obtained has been found before, by Lewis and 



Leach in a different but related context [^]. They found that a classical Hamiltonian with 
explicit time dependence has an invariant quadratic in the momentum if and only if the 



potential is of the form ( |4.19| ). (Their notation is related to ours via the substitution 
p(t) = 1/ f(t).) Lewis and Leach examined a purely classical problem, however; our inves- 
tigations are explicitly quantum mechanical. In this section, we show how to generalize f|| 
to the corresponding quantum mechanical problem. Remarkably, we find that a quantum 
mechanical Hamiltonian will have an invariant quadratic in the momentum in exactly the 
same cases that the classical Hamiltonian will, and that these are the potentials we have 
identified by the condition of separability. 

Any operator in quantum mechanics satisfies Heisenberg's equation 

Tt = m + n lH ' I] - <"> 

We are interested in the case for which there is an invariant I which is quadratic in the 
momentum. Because of the identity 

\p,A(x)) = -ih?^ , (5.2) 

we can always write / in the form 

/ = f 2 (x,t)p 2 + f 1 (x,t)p + f (x,t) . (5.3) 

Then one can easily derive that the functions fo(x,t), f\(x,t), and / 2 (x, t) satisfy the fol- 
lowing differential equations: 



dt 2m dx 2 m dx J 1 dx 

dt 2m dx 2 dx 2 dx 

Notice that in the limit H — > we recover the classical equations of Lewis and Leach 
as we expect. The above system of differential equations for the quantum case can be 
solved in the same manner as the classical equations. Remarkably, the solution to the 
system of quantum mechanical equations is virtually identical to the solution to the classical 
equations! In particular, we see that f 2 is independent of x, and that if we define fo(x,t) = 
f (x,t) — ihmf 2 (t)/2, then f , f%, and f 2 satisfy exactly the classical Lewis-Leach equations. 
We can use this to find easily that the form of all possible time-dependent potentials that 
satisfy the quantum Lewis-Leach equations is given by ( [4.19| ), and the associated integral of 
motion is 



1 L ( P 
2 



m 



pi a) — p(x — a) 



with p{t) = l/f(t), where f(t) is the function from the coordinate transformation ( [4.1| ). 
Note that this invariant is exactly the pseudo-Hamiltonian defined earlier in this paper. 

Of course, the existence of an integral of motion does not guarantee exact solubility; 
rather, it appears as the analogue of separability. Exact solubility will occur only for V(y) 
of certain functional forms. The shape invariance condition we used earlier in this paper, 
for example, does guarantee exact solubility of the corresponding time-dependent model. 



VI. PROPAGATORS 

The propagators for the potentials ( |4.19| ) can be given, of course, in terms of either the 
x or the y variables. The respective expressions are 

K(x, x'- 1, f) = £ * fc (x, t) %(xT, t') (6.1) 

k 
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and 



h 



(6.2) 



are obviously related by 



/ /x ^ fx — a(t) x' — a(t') , 
K(x,x';t-t') = K I 7-7^ — , -±J. ; t-lf 



(6.3) 



p(t) ' P (0) 

since y = (x — a(t))/p(t). (We find it useful to work in terms of p(t) = 1/ fit) in this section.) 
Using the wavefunction transformation ( |4.2j ) and substituting in the propagator expression 
JT2D, we find 

K(y,y';t,t') = (p(t)p(t')y 1/2 e*™ e*'™ £ S*(y,*) H^f) 



(p(^)p(i': 



-1/2 



e *(y,t) e 4>*(y',t') 



X 



Enfe/)ni!/')^ £t(TM • 



(6.4) 



Noting that for the ordinary Schrodinger equation ( 4.14 ) which arises for t), the prop- 
agator is 



K (y,y';t,t') = £ Y fe (y)nV) e"^^ 



(6.5) 



we see that 



K(y,y';t,t') = (p(t)p(t')) 



>\\~ 1/2 „*(V.*) „**(2/',t') 



^o(2/,2/';r(t)-r(0) 



(6.6) 



where r(£) is as defined in ( |4.17| ). Returning to the original variable x by using equation 
( p.3|) , and the results for $>(y,t) contained in (|4.6| )- (|4.8| ), we find, after some calculation, 
that 



K(x,x';t,t') = (p(t)p(i')) V2 exp (-^ J dth{t)\ exp 



isjC*(* (t >-$K 



x exp 
x exp 
xK \ 



' 1 JHl ( _ p^1 x '2 
2h \ P (t) P (t>) , 

t [ a{t) ~w) a( ) ) x ~ ir )_ W) a{ ] ) x . 

x — a(t) x' — a(t') 



pit) ' pit') 



7T^;r(t)-r(0 



(6.7) 
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Despite its length, this result is relatively elementary in nature. One should also be able to 
derive it in a straightforward manner in the path integral formulation J!]], although we have 
not attempted to do so. One can easily check that this propagator gives the correct result 
in the special cases a(t) ^ 0, pit) = 1 and a(t) = 0, pit) = ^rt' 2 + 2pt + q, by comparison 
with and 0, where these special cases are considered. 

We now give several examples based on the general formula ( |6.7|) . We will identify 
models in which our method in combination with other techniques makes the calculation of 
propagators relatively easy. In some cases we will actually display the propagators, while in 
other cases we will give the reader the necessary ingredients, but, in the interests of space 
and clarity, not present the full form of the propagator. 

As a preliminary, we rewrite the potentials ( |4.19| ) in the form 



V(x,t) = -^V^^y 1 ) -F(t)x\mu\t)x\ (6.8) 
where the functions u(t) and F(t) are solutions to 

pu 2 {t)p = , (6.9) 

au 2 {t)a = —F{t) . (6.10) 

m 

We will use u and F in lieu of a and p as convenient. 

The first example we consider is the free particle of time- independent quantum mechanics, 
i.e. V(y) = 0, which is of course an exactly soluble model. This implies that the forced 
oscillator with a time-dependent frequency 

V(x,t) = - F{t)x^mu 2 {t)x 2 , (6.11) 

is also an exactly soluble model. 

We can arrive at the same model by starting with a simple harmonic oscillator V(y) = 
muj 2 y 2 /2. Then by taking h(t) = — muj 2 a 2 /2p 4 we find exactly the potential Q6.ll ) with one 
minor modification. Now, p(t) is not a solution of (|6.9| ), but rather of 

pu 2 (t)p = 4- (6-12) 
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Equation (|6.10|) remains unchanged. This model is still separable and exactly soluble, even 
though p satisfies ( |6. 12 ) rather than ( |6.9| ), as the effect of the extra term in ( |6.12j ) can be 
absorbed into a re-definition of V{y) Q. 

It is known that the solutions of equations ( |6.9| ) and ( |6.12| ) are related. In partic- 
ular, if po(t) is a solution of equation fl6.12| ), then p\{t) = p (t) sinc^oTo (t) and p 2 (t) = 
Po (t) coscuoTo (i) , where r (t) = / * ds l/pg(s), are solutions of equation (|6.9| ) flOfl . Also, a 
particular solution for equation (|6.10| ) can be found using the theory of Green's functions: 

F(s) 



a(t) 



Po(*)Po(s) sin[c«; (ro(t) - r (s))] 



t 

ds 

to muj 

ds [p 1 (t)p 2 (s) - p 1 (s)p 2 {t) 

to muj 



(6.13) 



Alternatively, the above expression can be verified by direct substitution in equation fl6.10|) . 

Thus, given po(t), we can use either the propagator of the free particle or the propagator 
of the simple harmonic oscillator with constant frequency to obtain the propagator of the 
potential ( |6.11|) . Since the resulting expression is rather long and complicated, we will not 
write down the result here, as its evaluation is entirely straightforward. 

As a by-product, notice that the propagator of the time-independent oscillator of quan- 



tum mechanics can be derived by the one of the free particle as had been noticed in ||13|| . In 
this case, F(t) = 0, u(t) = ojq = constant, and it is enough to take a(t) = 0, pit) = cosujot; 
thus r(t) = (l/u;o)tan(c<;ot). Using the propagator of the free particle 

im (y — y') 2 



K (y,y';t,t') 



we find the well-known expression 



2irih(t' - t) 



exp 



2h t-t' 



(6.14) 



K(x,x';t,t') 



rn 



j imuiQ 

2TrihsmuJo(t' — t) ^ \ 2h 



(x + x )cota;o(£ _ 



2xx' 



cosu;o(i — t') 



(6.15) 



We now consider a related model, the free particle in 

0, if < y < L , 
+oo, if y < or y > L 



V(y) 



(6.16) 
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V{x,t) 



(6.17) 



The exact solubility of this model implies that the infinite well with moving boundaries 
—F(t) x + \ mu 2 {t) x 2 , if a{t) < x < p(t) L + a(t) , 
+00, if x < a(t), x > p(t) Lq + a(t) , 

is also exactly soluble. The propagator of the time-independent infinite well can be writ- 
ten down effortlessly using the energy eigenvalues E n = n 2 ir 2 h 2 /2mL/^ and eigenfunctions 



Y n(y) = J2/L sm(niry/L ), giving 



1 1 i\ 2 / i n 2 ir 2 h 2 \ mry nay' 

K {y,y,t,t) = — }^ exp - - 2 (t - t ) sin— sin— 



(6.18) 



This can be written in terms of the Jacobi ^-function, 

+00 

6 3 {u\z) = 12 Y, e mVz cos(2mx) 

giving 

K {y,y';t,t') 



(6.19) 



n=l 



2L n 



0, 



'n{y-y') 



2L 



irh(t - t'Y 



7T 



(y + y') 



2mLl 



2L 



2mLl 



(6.20) 



It is now an easy exercise again to write down the propagator for the infinite well with 
moving boundaries in (2, t)-space. 

Using the theory of supersymmetric quantum mechanics, Das and Huang fill have de- 
veloped a very elegant way of calculating the propagators of the time-independent exactly 
soluble potentials of quantum mechanics. For example, using this method one can obtain 
the exact propagator for 



V(y) = «L 2(2+1) 

2m y z 



n = 1,2,3,... 



(6.21) 



which is fTT| @ |3 



Our results, then, show that not only is the time-dependent potential 

rr/ x t^, x 1 9/ x 9 ^ 2 n(n + l) 

V(x,t) = - F(t)x-mu 2 (t)x 2 -— ^ — ^ , 

v ' ; w 2 w 2m (x - a(t)) 2 



(6.22) 



(6.23) 
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exactly soluble, but also that its propagator can easily be obtained from the propagator 
( |6.22j ) and our fundamental propagator formula (|6[ 



As final example, we note that in ||11|| , the propagator is found for the exactly soluble 
potential 



2fe2 



V{y) 



a 2 h 



rn 



^cosh (ay) 2 / 



(6.24) 



Putting the results of this paper together with those of ||11|| , then, we can easily find the 
propagator for the potential 



2fc2 



V(x,t) 



a 2 h 



cosh 



x — a(t) 



F(t)x\muo 2 (t) x 2 



(6.25) 



m p 2 (t) \ p(t) J 

The final expression for the propagator is lengthy, but simple to obtain. 

Needless to say, one can in the same way find the propagators for the time-dependent 



generalizations ( |4.19| ) of the other known exactly soluble time-independent potentials of 
quantum mechanics. The examples above have been selected as a sampler of typical results 
and techniques. 



VII. CONCLUSION 

In this paper, we have explored the applicability of separation of variables to quantum 
mechanical systems with explicitly time-dependent potentials. We have presented this work 
in one spatial dimension, but clearly the same methods can be applied in higher dimensions 
as well. We have found that only a very limited set of changes of variables and transforma- 
tions of the wavefunctions can produce separability. The uniqueness we find is remarkable, 
and it indicates that the underlying principles that allow separation of variables are quite 
restrictive. 

Using as our central technique separation of variables, we have also constructed exactly 
soluble time-dependent generalizations of any exact soluble time-independent model. In 
fact, our method reduces the identification of exactly soluble time- dependent models to the 
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identification of exactly soluble time-independent models. This allows us to make contact 
with the very powerful insights and techniques of shape invariance, and promises a deeper 
understanding of exact solubility in the time-dependent case. 

Our work also makes significant contact with other work on time-dependent Hamiltoni- 
ans. For example, we have used separation of variables to find the exactly soluble models 
discussed by f|. We have also worked out the quantum generalization of Q], and shown 
that the existence of an invariant quadratic in the momentum occurs for exactly the same 
time-dependent Hamiltonians as separability does. 

Of course, we cannot claim to have identified or classified all exactly soluble time- 
dependent potentials. However, the unique time-dependent generalizations produced by 
any useful changes of variables suggests that there is indeed a very intricate structure un- 
derlying the exact solubility of these models. To identify this structure — ideally, to do 
so directly in ordinary position space — would be a significant contribution to our under- 
standing of exact solubility and to a possible classification of exactly soluble time-dependent 
models. 

Furthermore, the results of this paper coincide at various points with the results of the 
references cited on exactly soluble time-dependent quantum mechanics 0-0], despite the 
differing approaches of all these papers (including ours). This suggests that separation of 
variables is a sound guide to finding exactly soluble time-dependent potentials. Given the 
variety of techniques that have been applied to time-dependent quantum mechanics, all of 
them of course producing mutually consistent results, there is probably much progress to be 
made by using these various approaches in tandem. 
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